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p_^- Abstract 
(D : 

^ . We construct a three flavor chiral Lagrangian of pseudoscalars and vectors witli special 

empliasis on tlie symmetry breaking terms. Comparing tree level two and tliree point 
I functions with experiment allows us to first, fix the parameters of the model (including 

^ ■ the light quark mass ratios) and second, to predict m{K*~^) — m{K*°), r{K* —>■ Kn) 

and r(0 — s> KK). The last mentioned quantities come out reasonably well, in contrast 
to an "ordinary" SU{?)) treatment. For this purpose we need "second order" symmetry 
breakers involving the vector fields analogous to those needed for the chiral perturba- 
tion theory program with only pseudoscalars. An improved description of the rj — r]' 
system is also given. We then use the soliton sector of this improved chiral Lagrangian 
to investigate some aspects of baryon physics which are especially sensitive to symmetry 
breaking. For this purpose a fairly elaborate "cranking" techinque is employed in connec- 
tion with the collective Hamiltonian. In addition to the "strong" baryon mass spectrum a 
careful investigation is made of the non-electromagnetic part of the neutron-proton mass 
difference. This work is needed to improve our previous estimates concerning the two 
component approach to the "proton spin" puzzle. We find that both the "matter" and 
"glue" contributions are small but they do tend to cancel each other. 
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1. Introduction 

In this paper we first discuss a three flavor chiral Lagrangian of pseudoscalars and 
vectors with special emphasis on symmetry breaking. The sohton sector is then used to 
calculate the non-electromagnetic part of the neutron-proton mass difference as well as 
the nucleon matrix element of the axial singlet current. Our motivation is to update some 
earlier papers [1,2,3] (which should be consulted for adequate references to background 
work). However, this updating has turned out to be non-trivial and to yield some results 
which may be of general interest. 

The proton matrix element of the singlet axial vector current has recently attracted a 
lot of attention because the EMC experiment [4] indicates that it is close to zero at = 0. 
This is called the "proton spin puzzle" [5,6] since the matrix element should equal unity 
(twice the proton spin) in the naive non-relativistic quark model. Amusingly, the simplest 
Skyrme model of pseudoscalars only does predict zero [7] . One might therefore take the 
point of view that the simple Skyrme model be considered a kind of first approximation 
to low energy dynamics. We have argued [1] that more is needed because additional 
"short distance" information is required to adequately explain the neutron-proton mass 
difference. In fact, the calculations of the n — p mass difference and of the axial proton 
matrix element are closely linked. Both require (at the dominant two fiavor level) the 
"excitation" of the r) meson field; the vanishing of one quantity implies the vanishing of the 
other. With vector mesons supplying the short distance component it is actually possible 
to explain [1] the neutron-proton mass difference and obtain [2,8] a suitably small value 
for the axial singlet matrix element. Another aspect of the proton spin puzzle which this 
kind of model might illuminate is the so-called "two component decomposition". This 
is an attempt to make the small value of the axial singlet matrix element plausible from 
the QCD parton model point of view. The original idea [9] of decomposing the axial 
singlet current into "matter" and "glue" parts has been criticized [5,10] as not being 
gauge invariant. However, a way to overcome this objection - by looking instead at a 
two component Goldberger-Trieman type of relation - has been suggested by Shore and 
Veneziano [11]. When this mechanism is implemented [3] in the meson Lagrangian it turns 
out to result, at the soliton level, in an extra contribution to the n — p mass difference 
proportional to the "glue" component of the axial singlet matrix element. A numerical 
estimate suggested [3] that the glue component, as well as the matter component, was 
small. 

Now, most of the above work consisted of calculations for the two flavor part but 
merely estimates for the effects of including the third flavor. In this paper we discuss 
the calculations using the full three flavor model throughout. It should be remembered 
that the n — p mass difference is driven by the symmetry breaking terms in the meson 
Lagrangian so it is especially important to treat these carefully. The same care is war- 
ranted for treatment of symmetry breaking in the collective Hamiltonian which describes 
the soliton sector of the theory. 

The flrst step of studying symmetry breaking in the meson Lagrangian is discussed in 
sections 2, 3 and 4. We want to have sufficient symmetry breaking terms to be able to 
explain at least the pseudoscalar decay constants and masses as well as the mass spectrum 
of the vectors and the 1~ — > 0~ 0^ decay amplitudes. How to do this for the pseudoscalars 
is well known. In the chiral perturbation theory (CPT) program [12] a list of symmetry 
breakers depending on the mass matrix in the underlying theory of QCD gives sufficient 
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flexibility to fit the 0~ decay constants and masses (apart from the 1] — rj' system for 
which a special treatment based on the U{1)a anomaly is required). There are also some 
loop diagrams resulting in "chiral logs" but these are numerically small for the indicated 
choice of scale . On the other hand, it has been generally felt that no such elaboration is 
required for understanding symmetry breaking in the vector multiplet. The vectors are 
typically considered [13] the most "normal" multiplet. Here we emphasize that this is not 
really the case. Both the isospin splittings and 1^ — > 0" 0^ decay amplitudes show large 
deviations from the simple SU (3) predictions. We find that these can be understood if 
suitable vector symmetry breaking terms which are the analogs of the dominant ones for 
the pseudoscalars are included. What it boils down to is the use of the symmetry breakers 
which satisfy Okubo's [14] original form of the quark line rule [14,15], except for the 0~ 
isosinglet channel. An improved fit for that channel is also discussed here in section 4. In 
section 7 we make some further remarks related to developing a larger analog of the CPT 
program which would also include vectors and presumably describe low energy hadron 
physics up to around 1 GeV. 

Returning to the main track of this paper, we describe, in section 5, the calcula- 
tion of the n — p mass difference using the improved meson Lagrangian. A number of 
improvements, which have been described in detail elsewhere, are included for treating 
the SU (3) collective Hamiltonian of the soliton sector. In the first place, the collective 
Hamiltonian is diagonalized exactly using the Yabu-Ando technique [16]. In addition, 
"cranking" corrections are included to allow for "centrifugal" effects with the rotating 
Skyrmions. This means that while the classical soliton involves the n"', uOo and fields, 
once it starts rotating, fields for 'q,r]' , K, K,u, p'^, and A'* get "excited" [17-20] with 
amplitudes proportional to the rotational "angular velocities" . This yields improved re- 
sults for "strong" baryon mass differences, static properties, etc. We find that the model 
which has been thus improved over that of ref. 1 continues to give a reasonable descrip- 
tion of the non-electromagnetic part of the n — p mass difference. What changes is the 
percentage of this difference due to the excitation of the rj and/or rj' fields. Previously [1] 
it was around 70% but is now reduced to about 20% for reasons discussed in section 5. 
The application of this calculation to determining the "glue" part of the proton's axial 
singlet matrix element is given in section 6. We continue to find that this contribution is 
small. However, because there is a smaller r] contribution to the n — p mass difference, 
the accuracy of this statement is actually decreased somewhat. 

2. Terms of the Lagrangian 

Here we collect together, for convenience, the a. chiral invariant, b. flavor symmetry 
breaking but quark line rule conserving and c. quark line rule violating terms of the 
pseudoscalar— vector effective Lagrangian. Some discussion will be given of why we have 
chosen to include the terms listed, but not others. 

The dynamical degrees of freedom are the elements of the 3x3 matrix of the pseu- 
doscalar nonet, (j) the 3x3 matrix of the vector nonet p^. We need the unitary 
matrices 

U^e^^'t'/^-, e = t/'/' = e**/^-, (2.1) 

where F^^ is a bare pion decay constant. The vector mesons, which also transform non- 
linearly under chiral U{3) x U{3), are related to auxiliary linearly transforming "gauge 
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fields" Aj; and by [21] 

Ai = ^p,^^ + kd,e, < = eV/^e + k^fi'^e, (2.2) 

where ^ is a bare pcpcp coupling constant. 

Note that we have included only the lowest lying s-wave qq bound states of QCD as 
our dynamical fields. 

a. Chiral invariant terms: [21-23] 

These include the kinetic piece for the vectors: 

- \Tr[F^y{p)F^^{p)\, F^^{p) = d^p^ - d^p^ - ig\p^, p^], (2.3) 
and pseudoscalar kinetic plus vector interaction pieces: 



ml{l + k) 



Tr{AlAl + + ^^^^^^^TriA^UAp^) 



8k V M M ■ M ■ 4^ 

= - -m^Tr(p^p^)-^rr[p,(a^ee^ + a,e^0]-^(i + ^)rK5.e5.e^) 

- ^(l-k)Tr{ed,e^d,0, (2.4) 

where is a bare vector meson mass and A; is a dimensionless constant chosen as k = 
{m.^/ Ft^qY to insure correct normalization of the pseudoscalar kinetic term. Please notice 
that we are using the "X4 = iV Euclidean-type metric convention. 

There are also important interaction terms proportional to the Levi-Civita symbol, 
^livap- To write these compactly it is convenient to use a form notation with one- forms 
a — dUU"^ and A^ and to write the action terms (rather than the Lagrangian density 
terms) : 

TwziU) + jTr [ici{A^a^) + C2{dA^aA^ - A^adA^ + A^aA^a) 

+ C3{-2i{A^fa+ ^A^aA^a)] , (2.5) 

where ci, C2, C3 are constants discussed in [20-23] while TwziU) is the Wess-Zumino-Witten 
term [24]: 



Eqns (2.3), (2.4) and (2.5) comprise a minimal (though nevertheless complicated in 
detail) model for pseudoscalars and vectors [25]. One might imagine extending it by 
systematically adding terms with two more derivatives. This would be analogous to the 
usual chiral perturbation theory (CPT) program [12] with only pseudoscalars included. 
It should be remarked, however, that a good deal of the four derivative structure of the 
pseudoscalar- only effective Lagrangian can be obtained by "integrating out" the vectors 
of a minimal-type pseudoscalar-vector effective Lagrangian. 
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b. Symmetry breaking terms: 

In this "updating" of section II of ref . [1] we will include only quark-line rule conserving 

terms, i.e. those which can be written as a single trace in flavor space. 

The current quark mass terms in the fundamental QCD Lagrangian may be written 

as 



jC-mass = -rhqMq, (2.6) 

where q is the column vector of up, down and strange quark fields, m = (m„ + md)/2 and 
M. is a dimensionless, diagonal matrix which can be expanded as follows: 

M=y\3 + T + xS, (2.7) 

with A3 = diag{l, —1, 0), T — diag{l, 1, 0) and S = diag{0, 0, 1). x and y are the quark 
mass ratios 

x^^, y^-^{ . • (2.8) 



rh 2 \ rh 

It is also convenient to define [13] 

R= ^--^ ^i^. (2.9) 
rud - m„ 2y 

The significance of R is that it can, in principle, also be determined from the rather well- 
known masses of the baryon octet. We have introduced the redundant notations above 
because of the importance of M. in (2.7), to which all the symmetry breakers are taken 
to be proportional. We then write for the quark-line-rule conserving but flavor symmetry 
breaking terms in the effective Lagrangian: 

Tr{M[a'{ApA^ + A^U^ A";) + ^'{d^Ud^U^U + U^d^Ud^U^) 
+^'{F^^uUF^, + F^^U^F;:^) + d'iU + t/t - 2)] 

+X'^[MU^MU^ + MUMU -2M^]), (2.10) 

where F^^^ — d^A^'^ — iglA^'^ , A^'^] and a', 7', 5', A' are the symmetry breaking 
coefficients. The parameters a", (5" etc. introduced in ref. [1] may be related easily to 
a', (3' etc. via: 

X = a" /a' = /37/3' = etc., y = a/a' = ft/ (3' = etc. (2.11) 

Our motivations for including the terms shown in (2.10) are as follows. First, the 5' 
term is the standard one which splits the pseudoscalar masses. Next, the a' term is the one 
with the correct chiral properties which splits the vector masses. However, as pointed out 
in ref. [1], the a' term also gives objectionably large derivative- type symmetry breaking 
for the pseudoscalars. This is partially canceled by the pure pseudoscalar derivative 
symmetry breaker given by the j3' term. The 7' term yields a derivative-type symmetry 
breaker for the vectors too. Note that the (3' term is the #5 type in the CPT classification 
[12]. The A'^ term, which involves two powers of Al, is the #8 type. It is argued to be the 
same order as the /?' term. Furthermore, it is the only remaining manifestly quark-line- 
rule conserving symmetry breaker to "second order" . From a practical point of view the 
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A'^ term also enables us, as we shall see, to vary x in (2.8) while keeping the pseudoscalar 
masses and decay constants fixed at their measured values. In the CPT program, the 
determination of A'^ is subject to the Kaplan- Manohar ambiguity [26]; we shall discuss 
this point later. 

c. Quark-line-rule violating and symmetry breaking terms: 

One of the remarkable features of low energy dynamics is the success of the OZI or 
quark-line rule [14,15]. There is however a conspicuous exception associated with the 
interactions of the pseudoscalar singlet particles {rj, rj'...); the physics of this is denoted as 
the "?7(1) problem". In a rough way, this picture is reflected in the CPT flt of ref. [12]. 
There the OZI rule violating terms of #4 type (~ Tr{d^Ud^U'^) Tr{M{U + W)) ) and 
of 7^6 type (~ [Tr{Ai{U + t^^))]^ ) are both claimed to be negligible. On the other hand 
the ^7 type OZI rule violating symmetry breaker (~ [Tr{A4{U — f/^))]^ ) is claimed to 
be non-negligible. We shall see that the latter is naturally associated with pseudoscalar 
singlet particle interactions. Incidentally note that we have now accounted for all the 
symmetry breaking terms involving only pseudoscalars which appear in the second order 
CPT list. To these we have added the a' and 7' terms in (2.10), which involve the vectors. 

Here, we shall write down essentially just those OZI rule violating terms needed to 
fit the 7] and rj' masses and two photon decay widths in an effective chiral Lagrangian 
framework. In order that the U{1) anomaly equation for the axial singlet current, df^J^ = 

G, with G = d^K^, = ^^^e^f3^sTr{Apd^As - ^^A^A^As) in terms of the 
QCD gauge fields, be obeyed in the massless limit it is convenient to introduce G as an 
"auxiliary" field in the effective Lagrangian. No kinetic term is to be written for G so 
it gets eliminated by its equation of motion. This procedure might be considered as an 
effective way of integrating over instanton field configurations. We thus add the following 
terms to our effective Lagrangian: 

-G^ + -^G \n{det U / det U'') + nTr {aa)Tr{a^ + iG e' Tr[MiU -U'')] 
K 12 

+ iGC'Tr[M{AiUA^-Ap^Ai)], (2.12) 

where K,n,e' and C are parameters. The first two terms in (2.12) are the standard ones 
for mocking up the U{1) anomaly and giving mass to the rj' meson [27]. The third term, 
which gives us the freedom to adjust the rj' "decay constant" , and the fourth term which 
is a symmetry breaker of "instanton-induced" type, have been mentioned before [28]. The 
fifth term is an "instanton-induced" symmetry breaker involving vectors; since it leads 
only to four and higher point interactions we shall neglect it for our present purposes by 
setting (' = 0. Now G appears only in (2.12). It therefore satisfies the equation of motion 

G^1x-'-^Tr[M{U-U% (2.13) 

where the mathematical ( SU (3) singlet) rj' field is defined by rj' = ^F^^Xi X being gotten 
from the decomposition U = e^^U with detU =1. G in (2.13) should be substituted 
into (2.12). Among other things it gives a term of #7 type in the CPT list. We have no 
special reason, on the other hand, to include a term of #6 type. 

An additional term involving G will be considered later in connection with an appli- 
cation to the "proton spin" puzzle. 
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d. Remarks: 



The total effective action of pseudoscalars and vectors is taken to be 

Feff = J [ (2.3) + (2.4) + (2.10) + (2.12)]c^^a; + (2.5). (2.14) 

Symmetry breaking terms proportional to e^j.^^ are being neglected here for simplicity 
(they do not contribute to most of the processes which will be discussed in the present 
paper). Similarly, OZI rule violation for particles other than the pseudoscalars will mostly 
be neglected. 

3. Physical quantities at the OZI rule conserving level 

We will find the parameters of the effective Lagrangian by comparing the tree level 2 
and 3 point functions with experiment. There is a natural separation into the properties 
of the pseudoscalar singlets (which we discuss in the next section) and the other particles 
whose properties follow just from the OZI rule conserving terms. 

Compared to the discussion in Sec. 11 of rcf. [1], there is the difference that the A' 
term is now included and also the 7' term is not set to zero. The latter results in a very 
substantial improvement in the predictions for the properties of the K* mesons and gives 
a different fit for the fundamental mass ratios x and y of the light quarks. 

First, expanding out the kinetic terms shows we should renormalize the fields as (taking 
typical examples): 

= ^K* Pl3ix, = Z^iPlllJ- + P22ix)/ V^, (t>n = P33^, (3.1) 

where 

= [1 + ^(^ - m]"\ Zk^[1 + ^(1 + x)(^ - 

Z,^Z^^{1- , = [1 - Ai{l+x)]"\ 

= (1 - 8x7')'/'. (3.2) 

Next, the computation of the "Noether" currents of this model leads to the identifi- 
cation of the physical pion and kaon decay constants F^rp and Fkp as: 

F-Kp = Zt^Ft^, Fkp = ZkF-j^ . (3.3) 

Expanding out the pseudoscalar mass terms yields, in the isospin limit, 

ml = ^{8' + ^X% m^ = -^[5'(l + a;)+2A'2(l + a;)2]. (3.4) 

7rp Kp 

Note that, as pointed out in Sec. 2, if A'^ were not present the two equations in (3.4) would 

lead to a fixed value for x = ms/m. However it turns out to be very useful to explore the 
physical situation for various values of x. The vector meson masses in the isospin limit 
are: 

ml = ml = {ml - 2a') /Z% ml. = K - 2a'(l + x)]/^l:*, 
ml^{ml-Axa')/Zl . (3.5) 
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Next, let us turn to the mass splittings between members of the same iso-multiplet. 
We are not computing the photon-exchange contributions here so the prediction will be 
for the non- electromagnetic part of the iso-splittings. We then find [| 

[m{K°)~m{K+)Un-EM = -=;^^[2ml{% - 4(3') - 46' - 16X'\l + x)], (3.6) 

[m{K°*) - m{K+*)Un-EM = yi-4^'ml. + 2a')/mK^Z],., (3.7) 

[Mp^]non-EM = (2a' - 47'mJ) /m^Zj. (3.8) 
(The p°—uj transition mass, Mp^ is defined in terms of the effective term in the Lagrangian: 

Finally, consider the pcfxf) coupling constants. The pTrvr coupling constant, gp-j^^^ may 
be defined in terms of the p 27r width, r(p — * 2tt) as 

r(p ^ 27r) = (3.9) 

where g(p) is the momentum of the final pion in the p rest frame. In our model. 

For the other vector decays into two pseudoscalars we find 

r(ir^ ^Kti) _ 3 ( mp F^p Zp Y \q{K^)\' 



Tip ^ 2n) 4 {rriK* Fkp Zk* ) \q{p)\^ ' 

T{ct>^KK) (mp (F^pV ZpV \qm^ 



r(p^27r) \FkpJ Z^j |g(p)|' 



(3.11) 



where q{K*), for example, denotes the |)seu(iosca/ar momentum in the K* rest frame. 
Also, the small OZI rule violation for decays was neglected in the second part of (3.11), 
as well as in the mass formula (3.5). 

Now let us consider the experimental determination of the various quantities intro- 
duced. This is actually the heart of the matter and contains several very interesting 
qualitative features. In the first place this determination yields estimates of the funda- 
mental quark mass ratios x = ms/rh and y = {rriu—md) / {mu+rrid) in our framework. The 
dependent quantity R = {1 — x)/2y is also relevant since it can be separately determined 
from consideration of other particle multiplets (like the baryon octet). For comparison 
we mention the results obtained in earlier work in this model [1] in which 7' = A'^ = 
and wherein Mp^ was the only experimental isospin violating quantity used for the fit: 

x = 37, ?/ = -0.36, R = 50. (3.12) 



^Only 3.8 was used for fitting in rcf. [1]. The extra factor ^ makes negligible difference. Note that 
the A'^ limit of (3.6) quoted, but not used for fitting, should be replaced by the present formula. 
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Let us also compare with the determination of Gasser and Leutwyler[13]: 

x = 25.0 ±2.5, y = -0.28 ±0.03, i? = 43.5 ± 2.2 . (3.13) 

A recent alternate approach to the CPT program by Donoghue and Wyler [29] predicts 
instead 

y = -0.54 ±0.09. (3.14) 

Those authors do not have a similarly precise estimate of x but consider R = 32 ±5 and 
hence x ~ 36 to be more acceptable. In the present paper we find a "best fit" for 

a; = 31.5, y = -0.42, R = 36 . (3.15) 

This is closer to (3.14) than to (3.13). 

The main qualitative difference between our result (3.15) and (3.13) can be roughly 
understood in the following way. Gasser and Leutwyler consider the vector meson nonet 
to be the most "normal" one in the sense that old fashioned SU{3) relations can be best 
trusted. Their value of R, for example, is essentially obtained from the vectors. However 
we find very important SU (3) violations for the vector nonet system. They did note 
that the non-electromagnetic part of [m{K°*) — m{K~^*)], which is related to Mp^^ (their 
basic isospin violation input) by SU{3), could not be adequately explained but attributed 
this to uncertainties in interpreting the experimental results. Here, we find that the 
experimental value for this mass splitting can be fairly reasonably explained if one allows 
non-trivial wave function renormalization for the K* and 0^ fields by keeping 7' 7^ 0. 
This also allows us to dramatically improve the predictions for the widths r{K*)/r{p) 
and r(0)/r(p) in (3.11). It is amusing to note that when SU{3) was first proposed, the 
symmetry prediction for r{K*)/r{p) worked very well. But since then, the measured 
width of the p has increased from about 100 MeV to about 150 MeV. The relatively large 
wave function renormalization for the K* tends to restore the agreement between theory 
and experiment. 

It is very easy to fit the symmetry breaking parameters x, y, a', 7', 6', X'^ as well as 
the "symmetric" parameters g and to experiment if one temporarily holds x fixed and 
computes everything else. First, feeding the known values for F^^p, FKp,m.„ and ttik into 
(3.4) gives 6' and A'^. From Fxp/F-^p = Zu/Zt, we next find the quantity (^ — 4/3') from 
(3.2). Equation (3.6) then yields y while the formulas for irip and m^. in (3.5) together 
with (3.8) give us a', 7' and m^. Next the p width via (3.9) and (3.10) gives g and finally 
(3' is found from (% —AfS'). The fitted parameters as functions of x are displayed in Table 
3.1. ^ 

We note that the strength of the standard non-derivative pseudoscalar symmetry 
breaker, 5', does not change much with x. The most dramatic effect is the increasing 
importance of 7' which, from (3.2), is seen to lead to an important wave function renor- 
malization for the K* and vector masses. When x = 25, the values of x and y agree 
with the Gasser-Leutwyler values in (3.13). As x increases it is seen that the strengths \(3'\ 
and |A'| decrease substantially. For our fit in (3.15), a' has also decreased in magnitude 

^Input parameters: F^p = 0.132 GeV,fVp = 0.161 GeV,m^ = 0.137 GeV, = 0.892 GeV,TOp = 
0.768 GeV,mA = 0.497 GeV, (m^" - mK+)non-EM = 5.28 MeV, Mp^ = -2.65 MeV, r(p -> 2tt) = 
0.1491 GeV,r( A'* ^ Kir) = 0.0498 GeV, r(</> ^ KK) = 0.0037 GeV. 



9 



X 


y 




a y\jiti V J 


(3'{GeY'') 


i 


5'(GeV^) 


|A'|(GeV") 


9 


iO 


— U. iO 


oo.o 


— y.uo X iu 


o c:n NX 1 n- 
— z.oy X iu 


-4 


— i.UO X iU 


-3 


Q 0/1 \x 1 n- 
C).z4 X iU 


-5 


1 /] cr XX 1 n- 
i.40 X iU 


3 


o.oo 


on 


n 1 Q 
— U. iy 


/I S Q 


/I 71 V 1 n— 3 

— 4. / i X iU 


1 /I o vx 1 n- 
— i.4z X iU 


-4 


/I Qo ^ 1 n- 
4.0Z X iU 


4 


Q 70 XX 1 n- 


-5 


c\ c^A XX 1 n- 
y.D4 X iU 


-4 


70 
O. ( Z 




n 07 

— U./ ( 


/I 8 
40.0 


IT? v 1 n~3 
— Z.UO X iU 


— / .yi X iU 


-5 


1 on X/ 1 n— 
i.oy X iU 


3 


Q no XX 1 n- 

o.yz X iU 


-5 


c Kc XX 1 n- 
0.00 X iU 


-4 


O.OO 


OS 


— U.oo 


4U.0 


S 70 v 1 
— o. ( z X iU 


tJ.OU X ±u 


-5 




3 




-5 


^ 1 ^ X 1 n~ 

tj. ItJ A lU 


4 


Q QQ 
O.OO 




—0 38 


38 3 


—2 25 X 10"^ 


-4.05 X 10- 


-5 


2.05 X 10- 


3 


4.01 X 10- 


-5 


4.31 X 10- 


4 


3 91 


32 


-0.43 


36.1 


3.43 X 10-^ 


-2.89 X 10- 


-5 


2.26 X 10- 


3 


4.04 X 10- 


-5 


3.51 X 10- 


4 


3.93 


34 


-0.49 


33.9 


8.46 X 10-^ 


-1.88 X 10- 


-5 


2.45 X 10- 


3 


4.06 X 10" 


-5 


2.71 X 10" 


-4 


3.95 


36 


-0.55 


31.6 


1.29 X lO--'^ 


-9.98 X 10" 


-6 


2.62 X 10- 


3 


4.07 X 10- 


-5 


1.84 X 10" 


4 


3.97 


38 


-0.63 


29.3 


1.69 X 10-3 


-2.19 X 10" 


-6 


2.77 X 10" 


3 


4.09 X 10" 


-5 


4.71 X 10" 


-4 


3.99 



Table 3.1: Parameters as functions of x 



X 


{K°* — K^*)non-EM 


T{p)/T{K*) 




m{(f)) 


15 


2.04 MeV 


5.U 


124 


1.01 GeV 


20 


2.44 MeV 


4.97 


103 


1.02 GcV 


25 


3.05 McV 


4.45 


81 


1.04 GcV 


28 


3.55 McV 


4.12 


67 


1.07 GcV 


30 


3.97 MeV 


3.90 


57 


1.09 GeV 


32 


4.47 MeV 


3.67 


47 


1.13 GeV 


34 


5.09 MeV 


3.44 


38 


1.18 GeV 


36 


5.84 MeV 


3.21 


28 


1.26 GcV 


38 


6.78 MeV 


2.98 


18 


1.43 GeV 



Table 3.2: Predictions 



and has reversed sign. Clearly S' and 7' arc the dominant symmetry breaking parameters 
for our fit. The value of the "gauge" coupling constant g does not change much with x. 

We now have four predictions which are given in Table 3.2: [m{K°*)—m{K^*)]non-EM: 
r{K* K'k)/V{p 27r), and r(0 KK)/r{p 2n). These quantities are sensi- 
tive to the wave function rcnormalizations of the vector particles containing strange quarks 
{K* and (p). First consider the K°* — K^* mass difference. The photon exchange contri- 
bution has been estimated to be —0.7 MeV [30]. The experimental value is alternately 
given as 4.5±0.4 MeV (if one simply subtracts the two numbers in the Review of Particle 
Properties [31]) or as 6.7±1.2 MeV (if one considers just the "dedicated" experiments). 
So we should have, 

[m{Kn - m(i^+*)]__^M « { . (3.16) 

This should be equal (see the hmiting forms of (3.7) and (3.8) ) in the SU{2,) hmit to the 
negative of the non-electromagnetic piece of the p — uj transition mass M^^ {non — EM) , 
which in turn has been found from analysis of experiment by Gasser and Leutwyler [13] 
to be — 2.65±0.20 MeV. Clearly the SU{3) prediction is very bad. However, agreement 
can be obtained in the present model if Zk* 7^ 1 in (3.7). We see from table 3.2 that 
the range of x between 34 and 38 is suitable for this purpose. This is due to the increase 
of 7' with increasing x. Next, consider the ratio of widths r(p — > 27r)/r(K* Kn) 
which experimentally is 3.0. Eq (3.11) shows that this ratio is also sensitive to Z^*- We 
see from table 3.2 that, once again, exact agreement is obtained for x around 38. Thus 
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considering just the properties of the K* particle one would be tempted to choose x =38. 
This represents (see table 3.1) a value of y significantly larger than (3.13) in magnitude 
but in agreement with (3.14). Larger values of x and y (compared to (3.13) ) are also 
favored by consideration of the ratio T{p 27r)/r((/) KK) which is experimentally 
40.3. We see from table 3.2 that best agreement for this ratio is obtained for x around 
34. On the other hand is fit best with a smaller value of x (about 20). While it is 
true that the simplification of neglecting u — (f) mixing has been made above, this effect is 
small enough so the prediction should be reasonably accurate. A compromise "best 
fit" with X ~32 (see (3.15)) improves three predictions remarkably, does not distort the 

prediction too badly and will be seen in the next section to be needed to lead to a 
good description of the 77 — 77' system. 

It is of interest to note that the vector meson wave function renormalizations for 
the best fit a;=32 case are = Zp = 0.99, Z^* = 0.84 and Z^ = 0.65. Since Z^, in 
particular 0, represents a rather non-trivial correction one might in the future want to 
investigate other possible higher order terms involving vectors and loop diagrams. In any 
event, our analysis shows that such large symmetry breaking is required to understand 
the vector mesons which contain strange quarks. 

We remark here that the additional symmetry breaking term: 

fi'Tr{A^MAijM) (3.17) 

can help to fine tune our results. Because it is quadratic in M. it will provide a factor of 
for m'^{(f)) and negligible contributions to other vector masses. Similarly, it will mainly 
affect the (p^ — > KK decay amplitude. There will be an extra contribution of — 2x^/i'/2'| 

to rn^icj)) and an extra factor of 1 + ^V^a for the KK amplitude. The choice 
/i' = 2.4 X 10-5 GeV^ then enables us to fit both m2(0) and T{p)/T{(j)) in Table 3.2 while 
still keeping x = 31.5. Notice that (3.17) is the analog for the vectors of the A'^ term 
in (2.10) for the pseudoscalars; there are then three analogous quark- line-rule conserving 
but symmetry breaking terms for each multiplet. 

4- Physical quantities for the ij — ij' system. 

In the simplest Lagrangian which can mock up the U{1) anomaly (just the first two 
terms of (2.12) included) the rj' mass and the rj — rj' mixing can be satisfactorily fit. 
However the 1] mass comes out too low 0. Here we will show that this problem can be 
remedied if the second two terms of (2.12) as well as the OZI rule conserving symmetry 
breakers in (2.10) are all taken into account. 

Defining the 2-dimensional vector r] = ( ] where 77^ = (0ii + (p22)/V^, Vs = 033; 
we collect the quadratic terms in the effective Lagrangian as 

-^d,r]^Kd,r]-^r]^Pr]. (4.1) 



■^In the extraction of the parameters ci and C2 of (2.5) from experiment in [22] some (j) decays are in- 
volved. We note that in (4.5) of this reference we should nowwrite e/Z^ instead of |e|. This replacement 

should also be made in (4.8) so that the old formulas for ci and C2 remain unaltered. 
*This can be seen by referring to Fig 2 of [32] . 
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Here the "kinetic" matrix is 



4 .a' 16 n 

Ktt = l + ;pi(-^-4/3') + ^, 

■K y TT 

Kss = 1 + t;^(^-4/30 + t;^, 

Kts = X5T=-^, (4.2) 



7r 



while the "potential" matrix is 

1 8 , K 



3 9 

P55 = -^[SxS' + ^XKe' + ^ + 8Kxh''' + 32x''X'% 
Pfi 6 is 



TT 



= P5T=^|/^6'(l + x) + ^ + 8«:X6'^]. (4.3) 

The "kinetic" matrix is diagonalized by 

R-\9,)K R{e,) = k = {^^^' IJ, (4.4) 

where 

/cos^ sin^A 

^ ^ V — Sm e*,; COS Oil ^ ' 



and Xi,2 = \[{Ktt + -f^ss) T \J {Ktt - KssY + ^K^s\- The angle 6*1 is obtained from 

tan2^i = -^^^. (4.6) 

J^SS — J^TT 

The effective potential matrix P' is then 

p' = k-^'^R-\e^)PR{er)k-^'''. (4.7) 

This yields the physical rj and 77' squared masses as the eigenvalues of (4.7): 

m'ir^.v') = \[{Ptt + P'ss) T \l {P^t - P'ss? + ^P?s\- (4-8) 

The "bare" column vector is finally expanded in terms of the physical fields r^p = ( ^, ) as 



T] = Rie^)k-^/^Rie2)vp, (4.9) 

where 61 is obtained from (4.6) while 62 is determined from 

2P' 

tan 2^2 = ■ (4.10) 

12 



Note that the transformation (4.9) is, in general, non-orthogonal. For comparison, the 
mixing convention in the octet-singlet basis is 



,^,=m[;.) (4.11) 

and, to the extent that K may be approximated by the unit matrix, 

6^61 + 62- 54.7°. (4.12) 

The standard source of information, which we shall also employ here, about mixing 

in the 7] — rj' system is comparison of the 7r° — > 27, 77' — >■ 27 and 77 ^ 27 decays. It is 
well known that 71° 27 is well-described by "gauging" the first term of (2.5). Thus it 
is sufficient to consider ratios of the other rates to the 7r° rate. We find 

r(vr° ^ 27) . T{r] ^ 27) _ r{r)' ^27) 9 _ , _ , 



ml m^, 2 



: {a,y : {a^Y, (4.13) 



wherein 

5 / cos 9i cos 6*2 sin 9i sin 6*2 
77^ £^ £^ 



5 / cos 9i sin 82 cos 82 sin ^1 \ / cos 9i co 

^ 72 V k\" ^ kl'' )^\^ 





Experimentally, {a^Y = 3.98 ± 0.60 and (02)^ = 7.20 ± 1.40. 

We may now try fit our parameters to these experimental data. Adopting the com- 
promise fit discussed in Sec.3 with x — 31.5, we have three new unknown parameters: 
K, e' and n (see (2.12) ). On the other hand there are four experimental quantities to be 
fit: m^,m^/, and the rj — rj' mixing information contained in (o-i)^ and (02)^ of (4.13). 
Conducting a numerical scan of the k, e' and n parameter space we find a best fit for 

K = 0.105 GeV^ n = 1.40 x 10"^ GeV^ e' = 1.59 x 10"^ (4.15) 



With these parameters the physical quantities are 

m(r)) = 549 MeV, m(r]') = 959 MeV, (a^Y = 3.98, (02)^ = 7.19, (4.16) 

in agreement with experiment. Furthermore, the quantities describing the non-orthogonal 
7] — T]' matrix in (4.9) are 

^1 = 7.44°, 62 = 34.7°, kl''^ = 1.07, kl'"^ = 1.36. (4.17) 

This would correspond, via (4.12), to a usual rj—r)' mixing angle of about —12.6°. However, 
because K ^ 1 this comparison is just suggestive. 

We have found that it is difficult to achieve a fit for the experimental quantities in the 
77 — 77' system for x significantly different from the value used above. It is also amusing 
to note that the A'^ symmetry breaker in (2.10) plays an important role in enabling us 
to get a good fit in the first place. Without the X'^ term present the only reasonable 
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fit was found witli very large and negative e' wliicli would give the wrong sign for the 
neutron-proton mass difference, as computed in this model. 

5. The neutron proton mass difference 

In the preceding sections we discussed the effective meson action with emphasis on 
the symmetry breaking parts, especially those which are isospin non-invariant. In this 
section we will use this effective mesonic model to obtain information about properties of 
baryons. In mesonic theories like the one under consideration baryons emerge as soliton 
solutions. 

In order to explicitly test the isospin non-invariant terms we will evaluate the non- 
electromagnetic contribution to the neutron proton mass difference: 

A = (iVf neutron ~ ^proton ) non- EM • (^•-'-) 

Experimentally the neutron proton mass difference is [31] 1.29MeV. However the elec- 
tromagnetic interaction (photon exchange) also contributes to the neutron proton mass 
difference. Using experimental electromagnetic form factors this contribution is estimated 
as [13]: 

Aem = (-0.76 ± 0.30)Mey. 

The negative sign agrees with Aem being dominated by the Coulomb repulsion. Thus we 
may extract: 

A = (2.05 ± 0.30)MeV. (5.2) 

a. Description of the approach 

It has been demonstrated previously [1] that vector mesons are a necessary ingredient 
in an effective meson theory in order for it to explain the neutron proton mass difference 
reasonably well. In ref.[l] however, only the SU{2) x f/(l) version of the model was 
properly calculated while the SU{3) part was just roughly estimated. It is therefore 
highly desirable to extend the treatment to SU{3). 

In the first step the static soliton with unit baryon number is constructed^ The 
appropriate classical ansiitze are: 

Substituting (5.3) into the Lagrangian C yields the classical mass Md = — / d^rC. The 
isospin breaking terms do not actually contribute to Md. An analytic expression for Md 
may be found in refs.[17]. Extremizing M^i yields second order non-linear differential 
equations for F(r),G(r) and uj{r). The solutions to these equations are characterized by 
their topological charge which is identified with the baryon number. The solution of the 
baryon number one sector has been extensively discussed in the literature[17-20,22]. The 
static soliton (5.3) is invariant under "grand spin" G = J+I transformations but not under 
spin (J) and isospin (I) transformations separately. Therefore the static soliton carries 
neither good spin nor good isospin quantum numbers. In SU{2) the projection onto good 
quantum numbers is carried out by introducing time dependent collective coordinates for 

^More details on the notation of this section may be found in refs [17]. 
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the zero modes of the theory. In SU(2) these are the spatial and isospatial rotations. 
These two transformations are actually equivalent due to the above mentioned grand 
spin symmetry. Thus in SU{2) the spectrum contains only baryons with equal spin and 
isospin, e.g. N,A. 

Unfortunately the extension of soliton models to include strange fields is non-trivial 
due to the presence of SU (3) breaking, i.e. flavor transformations arc not real zero modes. 
However, we will consider SU (3) to still be an approximate symmetry and therefore shall 
introduce collective coordinates for the whole flavor group. This approach allows us 
to easily make contact with ordinary baryon phenomenology. Furthermore it has been 
demonstrated in ref. [17] that the collective approach describes static properties of the 
low-lying ^ and | baryons reasonably well. An alternative point of view that S'f/(3) 
symmetry is strongly broken leads to the bound state approach[33]. In this approach 
a kaonic bound state is constructed in the background field of the static soliton. The 
corresponding bound state energy eigenvalue determines the mass splittings of baryons 
with different hypercharge since they are characterized by different occupation numbers of 
this bound state. Numerically the two (somewhat different) approaches yield comparable 
results for the mass differences of baryons with different hypercharge and spin. 

The collective coordinates A{t) G SU{?)) are introduced by flavor-rotating (5.3): 

i{v,t)^A{t)i,U^)i,A\t) (5.4) 

and 

i=p,(r,i) = A{t) [Pl^^^^ ^*^)A\t). (5.5) 



Here we have also allowed for a more general field configuration parametrized by and 
X*. We will see shortly how these additional fields get excited. The time dependence 
of the collective coordinates is made most transparent by the introduction of "angular 
velocities", (a=l,...,8): 



>lti=lj;AA = » + (5.6) 



wherein the Aa, (a = 1, ..,8) denote the SUiZ) Gell-Mann matrices. 

Obviously the vector meson terms in (2.4) as well as the "anomalous" part of the 
action (2.5) contain expressions which are linear in the time derivative. These expressions 
provide source terms linear in the angular velocity for the fields which vanish classically. 
Thus additional fields are excited by the collective rotation. Clearly they are linear in the 
angular velocities. The most general ansatz for the non-strange vector mesons excited by 
the isospin rotation is[17, 19-20] 

Po = ^[ei(r)0 + 6(r)(f-J^)f]-T, u;, = ^e,,.,fi,.f,. (5.7) 

The pseudoscalar nonet contains components which are excited by the isospin rotation[20] 
as well as by rotation into the strange directions [19]. Parametrizing: 
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suitable ansatze are rjT = \ixij') + Xsl^))? • O., rjs = \ixij') ~ '^Xs^ij'))^ " ^ ^-iid K = 
W{r)r ■ tQk- For later calculations it is convenient to introduce also an ansatz for the 
glueball field G = g{r)r ■ Q instead of eliminating G from the beginning via (2.13). From 
parity and isospin covariance we also get the following ansatze for the K* isospinor fields: 

K; = ^^K: K* = -l^^[tE{r)u + ^eij,rjn]nK. (5.9) 

Substituting (5.3-5,7-9) into the action and expanding up to second order in the angular 
velocities as well as linear order in isospin breaking yields the collective Lagrangian: 

L = -Mel + ^3^^ Q2 + 1^2 ^7^^ _ ^ ^3^^ ^^^^^ + ^7^^ ^^^^^ 

-|7(1 - Dss) - |7s(l - ^Is) - br E-=i D^iDsi - ^ts EI=^ D^^D^^ 
-FaDss - A3 ELMiDsi + r'ss^ss) + as E-=i ^3^^ + P3 EI=, D^Ma, (5.10) 

where the Dij denote the matrices ^Tr{XiAXjA^) of the SU{3) adjoint representation. 
The moments of inertia and are functionals of the radial functions (^1, ^2, ^, X, Xs) 
and {W, S, E, D), respectively. Varying the moments of inertia with respect to these fields 
yields linear inhomogeneous second order differential equations with the classical fields 
F, G and uj as source terms. Analytic expressions for the moments of inertia as well as the 
symmetry breaking parameters 7, ai and j3i may be found in ref.[17]. Here we only wish 
to explain the mechanism which excites %, x& ^"^^ 9- The non-strange combination x + Xs 
is excited by the terms proportional to e^i,ai3 in (2.14) exactly as explained in refs.[l,20]. 
The glueball field and the strange part of the 77 fields are subsequently induced via the 
equation of motion (2.13). 

The additional symmetry breaking parameters 75, 7^ and 7^5 are solely due to the 
(A')^ term in (2.10): 

75 = -IT = -^(A')'(l - xfldrr^l - cos2F), 

^^{X'f{l-xfJdrr\l-cosF). (5.11) 

The isospin symmetry breaking parameters T^, A3, and fS^ arc needed for the evaluation 
of the neutron proton mass difference and we present the somewhat lengthy analytic 
expressions in an appendix. Here we just wish to mention that and A3 contain only 
classical fields while and Ps contain expressions linear in the excitations (5.7-8) as well. 

The Lagrangian (5.10) is quantized canonically by introducing SU{3) right generators, 
Ra {a — 1, .., 8) via: 



dL 



-(a^^a + ttl^Sa + tt3^3a), a=l,2,3 

= -i^ = { -{f3^^a + PiDsa + P^D^,), a=4,..,7 (5.12) 



a=^ 



We separate the isospin breaking part of the Hamiltonian: 

8 

H^-Y. RaQa -L = Hi=o + Hj=i (5.13) 



0=1 



The isospin symmetric part Hi=q may be diagonalized exactly by generalizing the ap- 
proach of Yabu and Ando[16]. This yields the energy formula: 

E,., = + - pj(J + 1) - 3^ + ^ess, (6.14) 
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where ess is the eigenvalue of 

+P'^sil - Dl,) + P^^T Eil ^8^^8. + P'lTS EL4 DsaDsa, (5.15) 

and C2 = Z]a=i -Ra denotes the quadratic Casimir operator of SU{3). ess is obtained, 
of course, numerically, using differential operator realizations of the Ra [34]. The isospin 
breaking part is obtained to be: 

3 

= 13^38 + A3 5]p3*^8* + ^38^88) 
1=1 

+ ^ E D^^{R^ + «ll^8.) + 4 E ^3a(i?a + A^Sa)- (5.16) 

In the evaluation of the Hamiltonian only terms linear in the isospin breaking have been 
retained. Finally we have the non-electromagnetic contribution to the neutron proton 
mass difference: 

A = {n\Hi=i\n) - {p\Hi=i\p) 

3 
i=l 

-2^j2{p\D3i{R, + aiDsi)\p)-2^j2{p\D,^{Ra + f3iDsa)\p), (5.17) 

i=l P a=4 

where \p) denotes the exact proton eigenstate of Hi=q. In writing (5.17) we have made 
use of the flavor transformation properties of the D-functions. 

Although it is obvious, we would like to stress that in this treatment the neutron 
proton mass splitting is not obtained as the difference of two large numbers. On the 
contrary, the leading operator in the isospin breaking(if/=i) has been extracted and its 
expectation value is identified with the mass difference. 

b. Numerical results 

We now present our numerical results for the neutron proton mass difference. To pro- 
ceed we have to fix the parameters of the anomalous part of the action, (2.5). Previously 
[22] it has been shown that 

h = 4{2ci-C2/~g-C3/4f) ^OA 
~gvv<t> = 4c2^1.9 (5.18) 

fit reasonably well the decay processes ^ Svr and uj —>■ pvr. h and gwcp are allowed 
to vary ^in the range h = —0.15, ..,0.7 and gvv<t> = 1-3, ..,2.2 subject to the constraint 
\gvv<i> — /i| ~ 1.5 due to uncertainties in the determination of the u — (p mixing angle. 
The third parameter could not be fixed in the meson sector; however, it was argued that 

^Hans Walliser has pointed out to us that the experimental measurements for the relevant uj and 
partial decay widths have changed from the 1986 to the 1992 "Review of Particle Properties" . Continuing 
to use (4.4) [see footnote 3 above], (4.7) and (4.9) of [22] now gives us \e\/Z^ — 0.059 ± 0.005 and central 
values gvvct> = 1-81, h ~ 0.38. These are substantially similar to the older values. 
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X 




As^MeV) 






25.0 


-15.32 


0.08 


-0.0046 


-0.0171 


28.0 


-15.60 


0.05 


-0.0049 


-0.0178 


30.0 


-16.08 


0.04 


-0.0052 


-0.0185 


31.5 


-16.45 


0.03 


-0.0055 


-0.0189 


34.0 


-16.89 


0.02 


-0.0060 


-0.0195 


38.0 


-17.89 


0.00 


-0.0070 


-0.0205 



Table 5.1: The isospin breaking parameters T^, A3, and f]^ as functions of x. The input 
parameters are according to table 3.1. For parameters of the anomalous sector we take 
h = OA,gvv<i> = 1.9, K = 1.0. 



X 




Mi{MeV) 


M2{MeV) 


Mn - Mp{MeV) 


25.0 


-0.017 


-174.1 


-31.0 


1.60 


28.0 


-0.020 


-180.0 


-31.2 


1.75 


30.0 


-0.021 


-183.5 


-31.0 


1.88 


31.5 


-0.023 


-185.1 


-30.6 


1.97 


34.0 


-0.025 


-188.6 


-29.7 


2.12 


38.0 


-0.028 


-193.0 


-27.3 


2.41 



Table 5.2: The relevant matrix elements Mi = ^ T,f=i{D3i{Ri + aiDsi))p and M2 = 
-0 J2a=4:{Dsa{Ra + PiD8a))p to evaluate the neutron proton mass difference as functions 
of X. Input parameters as in table 5.1 

K = C3/2gc2 ~ 1 from studies of baryon properties. We will adopt this value together 
with (5.18). 

In table 5.1 we present the results for the isospin breaking parameters T^, A3, and 

as functions of x. Obviously A3 is negligible and we will omit it from now on. 

In table 5.2 we display the relevant matrix elements and the prediction for the non- 
electromagnetic contribution to the neutron proton mass difference, A in (5.3). 

Note that the 0:3 term dominates and contributes more than 50% of the neutron pro- 
ton mass difference. This is expected since the term represents essentially the pure 
two flavor contribution. As previously discussed [1] the use of the exact wave functions 
of the SU (3) Hamiltonian rather than the "unperturbed" wave function drastically sup- 
presses the matrix elements of the r3 piece. Clearly, reasonable agreement with the result 
extracted from experimental data (5.2) is obtained for x > 30. Especially the prediction 
for the neutron mass difference for the favored value in the meson sector a; = 31.5 turns 
out to almost coincide with the central value of the experimental data. 

Based on the discussion in ref.[17] one might however object that for h — 0.4 and 
gvv<s> = 1-9, the SU{2) moment of inertia = 4.52GeV^~^ is somewhat too small, 
predicting too large mass splittings between the i"*" and |^ baryons. A small variation 
to h = 0.36, gvv'p = 1-88 and a; = 28 yields = 5.00Gel^~^. The corresponding 
mass differences are displayed in table 5.3 and are found to agree reasonably well with 
the experimental data. Actually the SU{3) mass differences come out shghtly worse 
than previously since the effect of a' is lowered. This is due to the fact that large, 
positive «i and f3i terms in (5.10) mitigate the effects of a large collective symmetry 
breaking parameter 7. In the previous fit (without the 7' type term) a' was negative and 
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A 


S 




A 


S* 




Q 


Fit 


168. 


263. 


404. 


327. 


470. 


617. 


766. 


(5.18) 


144. 


237. 


350. 


360. 


480. 


602. 


723. 


Expt. 


177. 


254. 


379. 


293. 


446. 


591. 


733. 



Table 5.3: Best fit to the mass differences of the | and | baryons with respect to the 
nucleon for the parameter set h = 0.36, gvv<i> = 1-88, n = 1.0 and x = 28 (see also table 
3.1). Also listed are the predicted mass differences for the central values for h and gvv>i> 
with K, = 1.0 and x = 31.5 (3.15). 

significantly larger in magnitude than in the present case (cf. table 3.1). Noting]] that ai 
is essentially proportional to {l—x)a' we may understand why the SU{3) mass differences 
are increased. For this parameter set we obtain: 

A = 1.77Me1/, (5.19) 

which is slightly lower than the central value 2.05MeV of the experimental data for A 
but still within the error bars. The dominant contribution to (5.19) stems from the as 
term: 0.93Me^. Previously[l] we had found that this contribution was almost entirely 
due to the 1] excitation. This is now somewhat modified since firstly in ref. [1] only 
the unphysical ?7t was present and its mass of 550Mel^ was assumed. Employing the 
treatment of section 4 in the soliton sector as well decreases the effect of the 1] fields 
somewhat because they are suppressed by a larger mass. Secondly the 7' term which 
was not present in the approach of ref.[l] contributes about 30% to a^. The strange 
excitations contribute about 0.54MeV to the mass difference. The contribution of the 
term {0.29MeV) is strongly suppressed since the matrix element of D^s is much smaller 
than its SU{3) symmetric value when the exact eigenstates of Hi=q are used. Comparing 
(5.19) with the results obtained for the central values (5.18) we recognize a small change 
for the worse with the baryon best fit parameters. However, the better agreement of A 
for the set (5.18) is merely due to the smaller non-strange moment of inertia so we would 
expect the baryon best fit set to be more reliable for predicting baryon observables. 

In the next section we will discuss a small extension of the model Lagrangian which 
enables us to fine tune the predicted neutron proton mass difference by adding a term 
which may be interpreted as providing the gluonic contribution to the proton spin. 

6. Two component description of the proton spin puzzle 

In this section we will describe the implications of the preceding considerations on the 
so called proton spin puzzle which actually refers to the surprisingly small matrix element 
of the axial singlet current between proton states. In a naive quark model this matrix 
element measures twice the proton spin. As indicated previously [3] the neutron proton 
mass difference may be used to determine the gluonic contribution to this matrix element. 
In ref. [3] only an estimate of this quantity was made. Here we will present the complete 
calculation since all the necessary tools have been provided in the preceding sections and 
we may therefore disentangle the two components of the axial singlet matrix element: 
matter and gluon contribution [11]. 

''Correct charge normalization requires that we only include the contributions of the classical fields 
(5.3) to ai and f3i. 
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The axial singlet current, may most easily be obtained from 

4 = -2^ = V3F^d^7j' + sJl (6.1) 
with X being defined after eqn. (2.13). It is important to note that the second part, 

J', = e^.,,Tr{|(^ + |)p.p,K-v^72^F.,y,-4i^2(72 + 273)Kl^,K.} (6.2) 

is not a pure gradient and has therefore a non- vanishing matrix element at zero momentum 
transfer in contrast to the first part. In eqn. (6.2) we have used -^V^ = P^i — ^"^^j as well 

as p^{v^) = U^^"^"^ d^U^^"^ ± dfj^U^^'^U^^'^'^ . may be considered as the short distance part 
of the axial singlet current since it is dominated by the vector mesons and is not present 
in models with pseudoscalars only. The dimensionless constant s has been introduced by 
the replacement in the Lagrangian: 

d.xJl sd.xJl (6.3) 

to allow for a deviation from the nonet form for the pseudoscalar fields in the Lagrangian. 
The nonet form was originally introduced to satisfy the OZI rule; however such an ap- 
proach seems somewhat debatable for the rj' . 

Defining axial form factors of the nucleon for the flavor I — u,d,s via 

A/PoPn 

<P{p')\qn,7,qi\P{p)>= 



u{p')[^,^,Hi{q') + ^j,Hi{q')]uip), with q^^p,- p'^. (6.4) 



the relevant quantity for the axial singlet current is 



H{q')=jZm{q% (6.5) 

1=1 



'2\ 



Obviously the first term in (6.1) only contributes to the induced form factor H{q 
Yl^=iHi{q^). In the original Skyrmc model (without vector mesons) = and therefore 
also H{Q) = 0, which has been considered as a success[7] of the Skyrme model since it 
nicely describes the results of the EMC experiment [4]. 

In order to introduce the two component mechanism [11] it is necessary [3] to allow 
the pseudoscalar gluon field, G to couple to via the chirally invariant expression: 

-d.GJl (6.6) 

wherein i is a new dimensionless parameter. This additional term changes the equation 
of motion for the glueball field (2.13) to 

G= " 



+ (6.7) 
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s\t 


0.0 


-0.5 


-1.0 


-1.5 


-2.0 


0.8 


1.70 


1.84 


1.98 


2.13 


2.25 


0.9 


1.74 


1.87 


2.02 


2.16 


2.27 


1.0 


1.77 


1.91 


2.05 


2.20 


2.29 


1.1 


1.80 


1.95 


2.09 


2.24 


2.32 


1.2 


1.85 


1.98 


2.13 


2.27 


2.35 



Table 6.1: The neutron proton mass difference A as a function of the parameters s and t. 



We furthermore have the Ua{^) anomaly equation which is equivalent to the equation of 
motion for rj' after eliminating the glueball field: 

In both equations (6.7) and (6.8) we have neglected small symmetry breaking effects 
originating from the e' term in (2.12). Taking matrix elements of these equations between 
nucleon states we may read off the corresponding coupling constants at zero momentum 
transfer: 

s-t 2rnjv 

90NN - (6.9) 

t — S K 

Next we may use the coupling constants to reexpress the axial singlet current form factor: 

" ^^^^^^'''^^ " 27!F^'''^^^ " "matter" + "glue" (6.10) 

The modifications (6.3) and (6.6) of course enter the evaluation of the non-strange moment 
of inertia a^, however we find, even for a large deviation of s — t from unity, no significant 
change of a^. This is due to the fact that these modifications only affect the 77 and glueball 
fields but leave the vector meson excitations ^1,^2 and $ essentially unaltered. Since the 
7] fields do not explicitly appear in the definition of the corresponding matrix element 
numerically turns out to be independent of t; we find for the best fit parameters of section 
5 (see table 5.3): 

H{0) = 0.29s. (6.11) 

However, we do find a significant dependence of the predicted value for the neutron 
proton mass difference, A on s and t since there the rj fields enter crucially. In table 6.1 
we display our numerical results for A using the "baryon" best fit parameters. We note 
that for |t| < 1.5 the mass difference depends linearly on t while for larger values of \t\, 
non-linear effects which enter via the equations of motion are significant. 

To proceed, we consider the precise value of H{0) as the parameter which determines 
the "fudge factor" s via eqn (6.11). Then fixing t from the neutron proton mass difference 
allows us to separate the matter and gluon contribution to H{0) in eqn (6.10). We present 
our results for various experimentally allowed values of H{0) in table 6.2 for the parameter 
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H{0) 


t 


"matter" 


"glue" 




Baryon Best fit 


0.0 


-1.86 ±0.90 


0.54 ±0.26 


-0.54^0.26 


0.1 


-1.64 ± 1.02 


0.57 ±0.30 


-0.47 T 0.30 


0.2 


-1.31 ± 1.04 


0.58 ±0.30 


-0.38^0.30 


0.3 


-0.96 ± 1.00 


0.58 ±0.29 


-0.28^0.29 




Set (3.15,5.18) 


0.0 


-1.09 ±0.94 


0.32 ±0.26 


-0.32^0.27 


0.1 


-0.74 ±0.93 


0.31 ±0.30 


-0.21^0.27 


0.2 


-0.38 ±0.96 


0.31 ± 0.30 


-0.11^0.28 


0.3 


-0.00 ±0.97 


0.30 ±0.29 


0.00^0.28 



Table 6.2: The "matter" and "glue" contribution to the axial singlet matrix element of 
the proton for two sets of parameters. For the best fit parameters see table 5.3. 

set used to fit the baryon mass difference as well as for the central values of h and gvvti), 
(5.18) together with (3.15). We think that the results of the "baryon" best fit should be 
more reliable than the results for the set which gives a best fit to the meson properties. 
This is because for both the spin and isospin mass splittings, a fine-tuning of the non- 
strange moment of inertia turns out to be crucial (see section 5). In any event, the 
difference between the two sets provides a measure of the "systematic" uncertainties in 
our calculation. Taking into account that for the set (3.15,5.18) the absolute values for 
the gluon contribution to the axial singlet matrix element are lower than for the baryon 
best fit we might consider the results for the latter set as an upper bound. We should 
furthermore remark that it would be surprising if s were to deviate too much from unity. 
Thus, in each case the range for H{0) between 0.2 and 0.3 would seem most reasonable. 

To sum up, we may conclude from table 6.2 and the discussion above that while there 
is a tendency for the "glue" part to cancel the "matter" part, both are most likely to 
be on the small side compared to unity. This agrees with the results of [3], though it 
must be remarked that the present, more detailed, calculation allows larger uncertainties 
in the "glue" and "matter" pieces. The reason for this is that, as pointed out in the 
previous section, the dependence of A on the rj fields is lessened due to the inclusion of 
the 7' term in (2.10). This translates into somewhat larger values of (s — t) than were 
previously estimated [3] as well as large uncertainties in t. Since the present approach 
"fives off' the deviation of the calculated A from its "experimental" value, its accuracy 
could be improved if a more accurate value of the photon exchange contribution, Aem 
were available. 

7. Remarks on the mesonic Lagrangian 

The effective Lagrangian discussed in sections 2—4 has a number of connections with 
physical quantities outside the Skyrmc model approach. We would like to make some, 
partly speculative, remarks on this here. 

a. Light Quark Masses 

Their values are of interest, for example, in constructing models of the quark mixing 
(Kobayashi-Maskawa) matrix of weak interactions. Conventionally, [13] the absolute val- 
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md 


rris 


2.3 


5.6 


125 


3.2 


7.9 


175 


4.2 


10.2 


225 



Table 7.1: Light quark masses. 



ues of the (current) quark masses are considered to be "running" masses evaluated at the 
scale of 1 GeV. Estimates of the kind we have made here yield just their ratios. Evalua- 
tion of "QCD sum rules" yields [13] an absolute value for ms(l GeV) in the 125-225 MeV 
range. Combining this with the "best fit" in (3.15) yields the sets (all in MeV) given in 
table 7.1. 

Note that ratios like fjuid/nis = 0.21 are independent of the absolute value of 
m,(lGeV). 

b. Confirmation of R 

It is interesting to try to confirm the value R=36 in (3.15) from consideration of other 
particle multiplets. The most precisely known masses are, of course, those of the ordinary 
octet baryons. First order perturbation theory (which yields the fairly well satisfied Gell- 
Mann Okubo mass formula) gives the prediction: 

R= = ^ 7.1 

n-p E- -E° ^ ^ 

wherein each particle symbol stands for its mass, which in turn is assumed to have been 
"corrected" by subtraction of the photon-exchange contribution. Numerically this reads 
R = 60.7 ± 9.1 = 46.2 ± 5.9 . The large deviation between the two central values 
suggests a need to go beyond first order perturbation theory. If one expands around 
iTLu = rrid = rris = the procedure [13] is not straightforward since one encounters non- 
analytic terms like m^/^ and mlnm. A possible way out is to expand around the point 
f^u = T^d = fTis = |(m„ + nid + nis) rather than zero. It was recently shown [35] that this 
method leads to a second order prediction for R: 

2VSmT + {n-p) + (S° - S") 

where m-r is the "corrected" A — S° transition mass. In principle mj- could be determined 
by a precision measurement of the difference between the pK^ —>■ Arj and nK —>■ Arj 
scattering amplitudes, for example. In the absence of such information we can use ( [7.2|) 
to construct a lower bound for R if we assume that quantities at second order should not 
deviate too drastically from their first order values. This yields [35] roughly 

i? > 38 ± 10 , (7.3) 

which is evidently consistent with R=36. 



C. Kaplan- Manohar ambiguity 

If A4{x) is a 3 X 3 matrix field transforming as (3, 3*) under SU (3)l x SU {3)r, it is easy 
to see that [7W(x)]"'^^^ det7Vl"''(x) transforms in the same way. This suggests investigating 
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the substitution 

M = M' + bM'-^detM' , (7.4) 

where 6 is a real parameter and = TWMs the diagonal matrix proportional to quark 
masses, given in (2.7). Using the characteristic equation for a 3 x 3 matrix it can be seen 
that the effect of the substitution ( [r.4|) on the non-derivative symmetry breakers of the 
second order CPT program [12] is simply 

6'Tr [M{U + U^)] + X'^Tr{MtJMU + MU^MU^) 
+Ce[Tr{MiU + U^))]^ + C7[Tr{MiU - U^))]^ 

= 6'Tr [M'{U + U^)] + (A'2 - —)Tr{M'UM'U + M'U^M'U^) 

h6' ~ - h6' - - 

+ (C6 + ^)[Tr{M'{U + U^))f + (C7 + ^)[Tr{M'{U - U^))f + higher ordei(.7.5) 

Here we have chosen to use the octet chiral field t/, with det f/ = 1, for simplicity. By 
"higher order" we mean the terms proportional to h which are generated from the A'^, 
Cg and C7 terms. These are third order in the CPT program and are to be thrown away 
when one is working at second order. The significance of ( [7.5|) is that, assuming the higher 
order terms to be negligible, the identical physics is obtained using either the parameter 
set {m„, mrf, mg, A'^, Cg, c^} or the set {m^, m'^, m'^, (A'^)', Cg, Cy}; these are related by 

niu = m'^ + b m'jm'g , 
111(1 = "r^d + b ^'yjnils ; 
rUs = m'g + h m'^m'^ , 

A'^ = (A'7-Y' 
h6^ 
4 

C7 = 4 + ^, (7.6) 



C6 - 4 + -r ' 



4 

for any h. 

Evidently there is a one parameter family of ambiguities. A lot of attention has been 
focussed on the possibility of obtaining mu = in some "frame" since that might provide 
a non-axionic solution of the "6'-problem" . Assuming that a solution for the parameters 
is specified for the primed frame in ( [7. 61 ), we can transform to mu = with the choice 
h = —m'^/{m'^m'g). The corresponding transformation for R in (2.9) is then 

R^B!/{l+n^Jm',). (7.7) 

With the best fit choice (3.15) for the primed frame we find that i? ^ 26 is needed for 
rriu = 0; according to ( [7.3| ) this would perhaps be unlikely. 

The present model is, of course, different from the usual CPT program since, among 
other things, we are including the vectors. In any event, as discussed in section 2c, the 
assumption that OZI rule violation is dominant for the pseudoscalar singlet channel leads 
to a type 7 but not type 6 symmetry breaker. This would uniquely fix the "frame" for 
the present purpose. 
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d. Extension of the meson Lagrangian 

It is naturally of interest to consider how the present mesonic Lagrangian can be 
modified to yield a systematic continuation of the CPT program to higer energies. The 
CPT program [12] itself is most reliable for energies up to several hundred MeV, which 
restricts us to using just the pion fields. Each successive order of approximation adds terms 
with two more derivatives or one more power of the mass matrix (formally introduced as 
counterterms to loops computed at the previous order). In practice, going beyond second 
order may be difficult. Perhaps surprisingly, the three flavor continuation, including the 
K^s and the rj, seems to work reasonably well. However, it obviously cannot be extended 
in energy past the vector mesons without including them. Since we would like to use 
the meson Lagrangian to study baryons-as-solitons in the 1 GeV region, it would seem 
reasonable to add just the vectors (together, perhaps, with pseudoscalar and scalar singlet 
glue fields in order to develop a direct connection with QCD) and stop there. This provides 
a "clean break" in the sense of including just the lowest lying s wave states of the quark 
model. To model QCD in this way to all energies would likely require us to include all 
the higher resonances, as suggested either by their need to produce high energy Regge 
behavior [36] or by their presence in the large Nc approximation [37]. 

Once we have agreed to limit the "quarkonium" states to the pseudoscalars and vectors 
it is necessary to list the chiral invariant interaction terms. One question is whether this 
will introduce double counting since it is known [38] that "integrating out" the vectors 
reproduces the coefficients of some of the pure pseudoscalar terms. We think that this 
should not be a problem; one should include a priori both the vector and pseudoscalar 
terms and suitably readjust the coefficients of the pseudoscalar terms. For the "second 
order" treatment, loops computed with the "first order" terms ( (2.4) and the a' and 6' 
terms in (2.10) ) should be included. For a start, however, it seems reasonable to neglect 
them (while loop diagrams are of great conceptual interest they are often numerically 
negligible for the indicated scale choice 0) . 

Here we have included those terms which we believe to be dominant. Especially, we 
have neglected OZI rule violating terms, except as discussed in section 2c. The fit to 
experiment is reasonable for the 2 and 3 point functions considered. Trying to fit more 
quantities will probably require fine tuning by adding additional terms and including 
chiral logs. Considering the tremendous amount of data in the region below 1 GeV, it is 
clear that the best way to proceed is along an incremental "evolutionary" path. 
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Appendix 

In this appendix we list the analytic expressions for the isospin breaking parameters 
in the collective Lagrangian (5.9). 

Substituting the classical fields into the symmetry breaking part of the Lagrangian 

^For example, the "chiral logs" are not qualitatively important in (10.11) of the first of refs. [12]. 
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(2.10) yields: 



^ [r^cosF(F'^ + uj'^) - 2G{G + 2)cosF + 4(1 + G - cosF)] 

-2/?'cosF(F'V2 + 2sin^F) + 25'r=^(l - cosF) 
-^cosF(2G"2 + ^(G + 2)2 - 



2 , ./2„2^ 

12 



+^^^±^r\l-cos2F) (A.l) 



and 



A3 = —/={i - x)yX'^ / cirr2(2sin2F + 3cosF - 3). (A.2) 



For the isospin breaking parameters which appear in those terms which are Unear in 
the angular velocities we find: 



j drds Ps = Any j dr(3s 



Sttt/ 
a' 

ds = — {a;r2[2-2ei-cosF(2 + ei + 6)]-4$(l + G-cosF) 

+(X + Xs)smF[^u;\' + \f'\^ - G{G + 2)] 
-pcosFF'(x' + x's) + 2(x + X8)sinF]} 

-2(3'[{x + X8)sinF(F'V2 + 2sin2F) - 2(r2cosFF'(x' + X^) + 2sinF(x + Xs))] 

y2„2^ 



-t^[{x + Xs)smF{2G" + ^{G + 2f- u;"r' 



+4{uj'^y + 2G'$') + 2cosF(a;'(ei + Qr' + ^^G{G + 2))] 



-2S'r^{x + X8)sinF - 2eV2^sinF - s'^^X'^r^x + X8)sin2F (A.3) 



a' F 

(3^ = /iy(sin- + sinF)[cuV-2(l + G-cosF)2 
2gi^ 2 

F.,. „ . F ,^ F,, . F 



+2cos-[(l + G - cosF)(D + 2iysin^(l + cos—)) - ujr'^{S - 1 + cos—)]} 

+— {r a;sm — 
g 2 

F F F 

-W{1 + G- cosF)[sin— (1 + 2cos— + cosF) + 2sinF(cos— + cosF)] 

^%\^^(—F\E + 2sin|^W^' - F'W) - smF{{D + 2iysin|^(l + cos|^))]} 



+2(— - /?'){W^(sin- + sinF)(F'V2 + 2sin2F) 

F F..,,„ . F 



-2cos^(l + co^'^—)\r^F'W' + 2cos^sinFW^]} 
z z z 
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F F 
-4(5V^sin— (1 + cos— 

+^{cos^[G'i2D' + EG) + ^{G + 2)^D - 2uj'r\S' + '^E)] 



l+smF)[2G" + ^ 
-4A'V^W^[sin2F + 2sm^ + sinF + 2xsm9r(l + cos^r)] (A.4) 



-W{sm^- + sinF)[2G"' + ^(G + 2)^ - a;'^^]} 

3F 



Of course, for the parts of the Lagrangian aheady investigated in ref.[l] the expressions 
for Fa and coincide. 

Footnotes 

1. Only 3.8 was used for fitting in ref. [1]. The extra factor Z^'^ makes neghgible 
difference. Note that the A'^ — > hmit of (3.6) quoted, but not used for fitting, 
should be replaced by the present formula. 

2. Input parameters: F^^p = 0.132 GeV,F^p = 0.161 GeV,m7r = 0.137 GeV, mi^-* = 
0.892 GeV, rup = 0.768 GeV, tuk = 0.497 GeV, {mKo-mK+)non-EM = 5.28MeV,M^ 
-2.65MeV,r(p ^ 27r) = 0.1491 GeV, r(X* ^ Kn) = 0.0498 GeV, r(0 ^ KK) = 
0.0037 GeV. 

3. In the extraction of the parameters Ci and C2 of (2.5) from experiment in [22] some 
decays are involved. We note that in (4.5) of this reference we should now write 
e/Z^ instead of |e|. This replacement should also be made in (4.8) so that the old 
brmulas for Ci and C2 remain unaltered. 

4. This can be seen by referring to Fig 2 of [32]. 

5. More details on the notation of this section may be found in refs [17]. 

6. Hans Walliser has pointed out to us that the experimental measurements for the 
relevant cu and 4> partial decay widths have changed from the 1986 to the 1992 
"Review of Particle Properties". Continuing to use (4.4) [see footnote 3 above], 
(4.7) and (4.9) of [22] now gives us \e\/Z^ = 0.059 ±0.005 and central values gvv4> = 
1.81, h — 0.38. These are substantially similar to the older values. 

7. Correct charge normalization requires that we only include the contributions of the 
classical fields (5.3) to cci and Pi. 

8. For example, the "chiral logs" are not qualitatively important in the first of refs. 
[12].). 
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